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ABSTRACT 

Here proposed a modification of the Press-Schechter theory allowing for the presence of 
a host density perturbation - host halo or void. The perturbation is accounted as statistical 
constraints in a form of linear functionals of the random overdensity field. Deviation of the 
background density within perturbation is interpreted in a pseudo-cosmological sense. Result- 
ing mass function of sub-haloes depends on the perturbation parameters: its mean overdensity, 
spatial scale, and spatial momenta of higher orders. Applications of the theory to superclus- 
ters, voids and bias problem are briefly observed. In its present form, the theory can describe 
the clustering properties of sub-haloes inside a non-virialized host only. Possible fix of this 
drawback is also discussed. 

Key words: large-scale structure of Universe - galaxies: luminosity function, mass function 



1 INTRODUCTION 

' The theory of Press and Schechter originated as a semi-analytical 
, approach to describing the evolution of the mass function of dark 
haloes. Today this theory, together with its modifications, perhaps 
the only one that, using some reservations, most closely matches 
, observations and numerical modeling for the widest mass range of 
' dark haloes, down to the resolution limit of numerical runs. The 
reservations of the theory concerned to its main provisions : 

- reinteipretation of halo merging as a random walk process al- 
lowed to complete the const ruction of the theo ry with a mass func- 
tion of Press and Schechter (iBondetal .1119911) : 

- iBenson et alj j2005l) showed, in the original theory the merger 
rate, or the merger kernel, in terms of Smoluchowski equation, is 
asymmetric function of mass which is flaw; the procedure to find a 
symmetric kernel proposed by the authors can solve this issue; 

- spherical collapse model used to get an overdensity threshold 
seems too rough; inste ad were adopted the model of ellipsoidal col- 
lapse ( lMo naco"l997a''^ and more general non-spherical collapse 
model {Tee & Shandarin 1998). 

Other reservations related to the environmental effects. Namely, to 
account to the prescripted large-scale distribution of haloes or the 
super-/sub-halo relation, w ere proposed mode ls using considera- 
tions about merger process dMo & Whitell 19961 on alternatives see 
references therein), and also heuristic models dPeacock & SmithI 
[2000.) . Q Modifications to the Press-Schechter theory proposed in 



* E-mail:kurbatov@inasan.ru 

^ Please excuse the author for the reference list is far fro m completeness. 
There are over 2000 citations of .Press & Schechteilfl974 in ADS, it's hard 
to analyse all of them. 



this paper allow for the presence of a host density perturbation - 
host halo or void, and give the mass function of sub-haloes depend- 
ing on the perturbation parameters. 

Below in 2nd Section we will describe a modification of the 
Press-Schechter theory. In 3rd Section we consider applications of 
the theory to superclusters and voids. Benefits and issues will be 
discussed in 4th Section. 



2 PRESS-SCHECHTER THEORY WITH STATISTICAL 
CONSTRAINTS 

2.1 Outline 

We assume that the basic model of the theory is canonical, i.e. 
cumulative distribution function for variance S is set by the ran- 
dom walk in overd ensi ty space with variance acting as a 'time', see 
iBond et alj ( Il99lh and lLacev & Coi3 h993l) : 



F{< S) = erfc 



.5, 



(1) 



where <5, is the threshold overdensity for the random process, a 
barrier for trajectory to pierce to be associated with a collapsed 
object. Corresponding differential probability distribution function 
(PDF) is 



fs 



dF{< S) 



5, 



exp 



Si 

'2S 



(2) 



dS 

ISheth & ToniienI ( Il999h suggested a correction for this formula 
which gives better fit to their A^-body simulations data but we will 
use the canonical variant for clarity. Mass PDF, /m, is obtained 
by differentiation F by mass through the mass-dependent variance. 
Both, threshold overdensity and variance depend on time via the 
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cosmological model, power spectrum of fluctuations, and model of 
collapsing lumps. Defining these dependencies we get an evolution 
of the mass spectrum of haloes. 

The value for threshold overdensity is commonly assumed to 
be 1.69; it is obtained as a critical overdensity for collapse by con- 
sidering the growth, turnaround and collapse of a uniform sph eri- 
cal overdense region JSardeen et al.lll986l:lLacev & Cold 19931). In 
general it depends on cosmological density parameters feke et all 
1 19961) but remains close to the conventional. Parameter 5, can be 
interpreted not only as collapse condition but also as a marker for 
mass of the 'nonlinear' structures. Since the variance is the mono- 
tonic function on mass, PDF reaches its maximum value at the 
unique mass which sense is the characteristic or nonlinear mass, 
m*. The extremum point for PDFs /g and fm satisfies equa- 
tion iSijn,) — 5^. On the other hand, in evolved systems the 
characteristic mass is th e mass of a most frequent occurrence. As 
IZel'dovich et al. I ( Il983h noted, the variance of the density grows 
rapidly at this mass and crosses unity, so we can reinterpet 5, 
as an overdensity in objects just evolved to the nonlinear ampli- 
tude. For this we assign threshold overdensity to a value marking 
the c haracteristic mass (qualitatively the same was mentioned by 
llain & Bertschinger C1994,) ). 

5* = ^/3 



1.73 



(3) 



This value is quite close to commonly used 1.69. Later we will 
assume Eq.[3]against the conventional. 

Connection of the variance to the mass of the collapsed struc- 
tures is implemented by filtering the local variance of density fluc- 
tuations with some kernel, which is considered isotropic usially. 
A space scale of the kernel, by-turn, is bound to the mass. Since 
fluctuations' amplitudes were small at early times, the binding is 
simply m — (47r/3) ilmPo^Jf , where R{ is the filtering scale, and 
flmpo is the mean matter density in the Universe. Denote / the 
random overdensity field. The filtered one in the real space is 



f{t,r,Rw) 



d\'Wir-r',Rw)f{t,r') 



(4) 



where W is the filter; the filtering scale _Rf is bound to _Rw with 
some filter-dependent relation, see Sect. 12.31 For the given random 
field the variance is 



S - 



{7-{7)y 

d^kd^k' W{k, Rw) W{k', Rw) K(t, k, k') 



(5) 



where W is the Fourier image of the filtering kernel; K{t, k.k') — 
{f{t, k) f{t, k')} is the correlation function for Fourier modes of 
unsmoothed field; the proportionality sign means free choose of a 
normalization constant of the Fourier transformation. The field is 
often assumed statistically uniform, isotropic, having zero mean 
and independing Fourier modes, i.e. Gaussian. This reduces the 
variance to depending only on isotropic power spectrum P in mo- 
mentum space: 



S(xJ d^kW^{k,Rw)P(t,k) , 



(6) 



where K{t,k,k') = S-D{k — k') P{t, k), i.e. modes are delta- 
correlated. With regard to general models, Eq.|5] if non-Gaussianity 
has form of the statistical constraints applied to the Gaussian field, 
then the correlation function can be obtained completely from 
power spectrum. Derivation of the correlation function in this case 
will be done in next Subsection. 

The power spectrum itself depends on the cosmology and is 



represented usually in a form of the primeordial power spectrum of 
the given spectral index n modified with transfer function T: 



P{k) oc fc"T^(fc) 



(7) 



As an example of the transfer funct ion can be mentioned the ap- 
proximation of lBardeen et al.l ( ll986h : 



T(fc) 



ln(l + 2.34q) 
2.34q 



X [l + 3.89g+(16.1g)^ 



(5.46g)^ + (6.74g)*] 



(8) 



where q = fc/F, and F is the shape parameter. Transfer func- 
tion can also be computed numerically as the result of the evo- 
lution of perturbations in an early times, e.g. by CMBFAST code 
(Seljak & Zaldarriaga 1996). The unconstrained variance calcu- 
lated for this power spectrum then must be normalized to the given 
value of cr| = S{R{ = 8ft"^Mpc). 

Calculations for evolution of the power spectrum is a difficult 
task even in the lowest orders of the perturbation series approach. 
For Gaussian field these calculations were completed just up to 1- 
loop corrections, i.e. to 2-n d order of accuracy in power spectrum 
jjain & Bertschingeilll994l) , or up to 3-nd orde r in expressions for 
filtered statistical momenta jScoccimarrol 19981) . In a non-Gaussian 
case the task became more complicated since the field is not zero- 
me an, and the same precision order requires much more integrals to 
get dCrocce & Scoccimarroll2006bllal) . In our theory we will restrict 
ourselves to a purely linear evolution, so the overdensity will be 
proportional to the linear growth factor D{t): 



f{t) = D{t)f^ , 
K{t) = D^{t)Ki^ 



P{t)^D\t)Pi^ , 
Sit) = D\t)SL , 



(9) 



where L index means values linearly evolved to the present time 
with unity growth factor. 

The linear evolution of perturbations settled on a large scale 
host overdensity (of both signs) can be represented in a linear per- 
turbation theory for the certain cosmology. Parameters of such a 
pseudo-cosmological model determined by the mean density of the 
host: an overdensed background corresponds to higher value of the 
matter density parameter in the pseudo-cosmology than in the Uni- 
verse, and vise versa. In the overdense host the perturbations will 
attain higher amplitudes than in global cosmology, in the under- 
dense one the last will be lower. Since D is the only factor provid- 
ing the time dependency in the class of the Press-Schechter the- 
ories, the bigger values of it will be in charge of the greater age 
of objects population. As the population evolves mainly by merg- 
ing of the objects, there is a natural connection between the host 
overdensity and the intensity of a merging process for sub-haloes 
of different mass. At this step we should obtain the bias factor for 
the halo number. 

The result of the theory should be the number density of 
objects of a given mass in a predefined host halo, and also sec- 
ondary values like cumulative mass of a population, etc. The num- 
ber density per unit mass interval is usually defined as drim/dm = 
rii-nPofm/m. In this definition is assumed the spatial homogenity 
of the distribution, i.e. there is no spatial structure inside the host 
perturbation. It might not be so in a virialized host halo, where the 
merger process and therefore the mass function can be affected by 
the shaip features of the distribution. This may lead to a mass seg- 
regation of sub-haloes in space and, hence, to redistribution of the 
mass spectrum of sub-haloes. We assume this feature can be mod- 
elled by a mass-dependent factor cjrn , so the general definition for 
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the number density is 



dm 



frr 



(10) 



Factor ujm obviously makes a contribution to the bias factor. Other 
distribution functions should be represented similarly. 

In this work the ACDM model is used with following param- 
eters: Qa = 0.7, ilm = 0.3, £78 = 0.9, spectral index for power 
spectrum is n = 1. The transfer function was computed using the 
CMBFAST code of lSeliak & Zaldarriagal (ll996V The units adopted 
are/i^^M© for mass, /i~^Mpc for length and (hHioo) ^^ for time. 



2.2 Constrained correlation function for modes 

We represent the host halo in terms of a functional constraints for a 
random field of the overdensity perturbations. In our task it's con- 
venient to make calculations using spherical modes decomposition. 
To obtain the correlation function for amplitudes of the spheri- 
cal modes we will use the approach of Hoffman & R ibak (1992), 
which was proposed by them for plane waves. 

The field /(r) = f{r, 9, (j>) can be transformed to spherical 
modes via the transformation 



fir) 



EE 



dkkji{kr)Yim{0,(t>)fin.{k) , (11) 



where ji is the radial Bessel's functions, Yim is the spherical func- 
tions with normalization dQ Y*^Yii^i — ^n'^mm' • Amplitude 
of spherical decomposition of the field, or image, is 



dr^kji{kr)(t dnYiUe,ct))fir) . (12) 



Transformation for radially symmetric field H{r) consists of only 
the isotropic modes: 



Hlm(k) = Sio5moH{k) , 

where Sa is the Kronecker's delta, and 



H{k)^V4^J- drr^kjo{kr)H{r) 



while 



1 



dkkjo{kr)H{k) 



(13) 



(14) 



(15) 



/47r V TT J„ 

Let us write constrains in a form of linear functionals fixing up 
a value of the convolved field at a point. Since the local extremum 
is the only constraint we use, we adopt the radially symmetric con- 
straining kernel, and place the point at origin. Using this conditions, 
the functional for a-th constraint can be written in form 



C,[/l= / d\//(„)(r)/(r) 



dfcH(,)(fc)/oo(fc) . (16) 



Constraints itself are fixed by assigning values to these functionals. 
Actually kernels can depend on a set of parameters characterizing 
a host halo, like a space scale or some momenta. We will not write 
them as argum ents for shortness til l next Subsec tion. 

Following iHoffman & RibakI jl99lL [l993) it can be showed 
that the pair correlation function for the constrained ensemble of 
modes is 



Klml'm'{k,k') = (^flm{k)fl>m'{k'^{Ca} ^ 

= ((flnik) - ^t'(fc)) (fl'n.'{k') - f^tik') 



(17) 



where the averaging performed over the unconstrained ensemble. 
Here 

fl^ik) = SioS^oQ'l CM] H(,a)(fc) (18) 

is the image of the mean constrained field but with non-fixed func- 
tionals C'a[f]', H(^)(fc) is the image of cross-coiTelation function 
between the field and a-th constraint: 

(T-) = {/ [/]) , H(,) (fc) = ^(,) (fc) P{k) ; (19) 

and Q^p is the inversion of the constraints correlation matrix: 



Qcp = (Cc [f] Cfi [/]) = j dk (fc) H^p) (fc) P{k) . (20) 
The unconstrained field is delta-correlated, which means 

{flr^{k)fvra'{k')) =5u'5rr,rr.'5'i\k-k')P{k) , (21) 

where is the one-dimensional Dirac's delta-function. Our goal 
is to find the variance of the constrained field, i.e. the full convolu- 
tion of a kind Eq.[5]for coiTelation function KimVm' (k, k') with a 
filter W . It's easy to show that the only part of last meaning is 

K{k, k') = 5^^^ {k~k') P{k) - Q-l H(„) (fc) S(^) (fc') . (22) 

The filtered variance can finally be expressed as 

S = Jdkdk'W{k,Rw)W{k',Rw)K{k,k'). (23) 

Note that unfiltered convolution implies using the image for 
three-dimensional Dirac's delta-function as a filter, (5^ ' (fc) — 
k/{2^^^n). In general, the filter can be written in the form 
W{k,Rw) = <5^^'(fc)m(fci?w), where iu(0) = 1, w(cxj) = 0, 
and ^5 juij ^1. The same applies to the constraining kernel. This 
is the result of the definition for unconstrained correlator, Eq. |2T] 
where the one-dimensional Dirac's delta-function is used. Let's set 
a single constraint whose kernel is the filtering kernel also, and the 
constraint's parameter is a space scale _Rh. The variance is then 



5* = 



J dkU{k,Rw,Rw) 



1 



W{k,RH) J dk'U{k',Rw,Ru) 



W{k,Rw) Jdk"Il{k",RH,Rii) 



(24) 



where n(A:,i?w,-RH) = W(fc, i?w) W(fc, i^n) P(fc). It is obvi- 
ous that the expression in braces tend to the values close to unity 
when Rw «C Rh or Rw 3> -Rh, and vanish if _Rvv ~ Ru- There- 
fore, we can expect, in general, that the abundance of the haloes 
vanishes if halo mass tends to the mass of the host perturbation, 
and came closer to its unconstrained value if halo mass a lot differs 
from the host mass. It is not surprising. 

The spatial halo correlation function is 



ar) = {f{0)f{r)\{C^}) 

= ^E / dk5^^\k)ji{kr)Y,r, 



(25) 



(26) 



X / dk'5^^\k')Ki^oo{k,k') 



For a spherically symmetric consraints it's image is 

C>) = 5<,="(fe)P(fc)-Q-^H(„)(fc) / dk'S^^\k')E^^^{k') 



(27) 

Time dependence of the variance and correlation function is 
determined by assuming linear growth, Eq. |9] i.e. by renormaliza- 
tion of the variance with the factor (see Subsection l2.4t . 
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2.3 Constraining and filtering kernels 

Let H{r ) be the kernel for the functional constraining peak (or dip) 
amplitude of a iield at origin. Besides amplitude, the shape can be 
constrained also via derivations of the kernel: 



{1st deriv) 



dH{r) 



H. 



' --(2nddcnv)\. , - 

and so on. Using reccurence relations for radial functions ji and 
their derivations it's easy to show that constraints for odd-order ra- 
dial derivations turn to zero because of symmetry. Images for the 
second and fourth derivations are expressed as 

-H'{2nd deriv) (fc) = ^ y H (k) , -H'(4th deriv) (fc) = y H (k) . (29) 

Filtering kernels generally accepted to use are k-sharp. Top- 
hat and Gaussian: 



(r) = 



d2ff(r) 



(28) 



(30) 
(31) 

(32) 
(33) 

(34) 
(35) 



- k-sharp 




WKs{r,Rw) ■- 


1 ji{r/Rw) 
2-kR%j r/Rw ' 


WKsiKRw) 


^5l^\k)e{l~kRw) ; 


- Top-hat 




WTHir,Rw) 




WTH{k,Rw) 




- Gaussian 




WG{r,Rw) = 







Wa{k,Rw) = 5^j,'>{k) e 

where ~5'i\k) = k/{2^''^ n) is the image of the three-dimensional 
Dirac's delta-function. Consider one of these kernels as the con- 
straining kernel which depends on only the single parameter, a 
space scale _Rh. Volume integral for any of these filters is normal- 
ized to unity, although it is not necessary for constraining kernels 
since any factor will be reduced in Eq.|22] However, this normaliza- 
tion allows for clear physical treatment of constraints. For example, 
the mass of the host perturbation with overdensity A'*^' (r) inside a 
sphere of radius r is 



M(<r-) = J p{r')e{l-r'/r) 



(36) 



(37) 



where the second term in braces is the overdensity averaged in the 
volume of the host. The host overdensity profile can be obtained 
using Eq. [TS] for the fixed values of the constraining functional 
Ca, then the averaged one will be 



(38) 



A^Q-'^Cc. / dfcM^TH(fc,r)i?(^)(/c,i?H)P(fc) 



Relation between the characteristic scale Rw and filtering ra- 
dius Rf, as well as between the scale Ru and actual constraining ra- 
dius _R h is deteiTnined by the choise of the kernel. As Bardeen et al.| 
( Il986l) noted, the characteristic scale for the Top-hat filter is not 
the actual space scale which is filtered out if we interested in the 
mass enclosed inside, but connected to it via Ri — 0.64i?w. For 
this ratio and specified power spectrum normalization, erg = 0.9, 



the nonlinear mass for unconstrained power spectrum is m* « 
3 X lO^^/i~^M0. This value is used to adjust characteristic scales 
for k-shaip and Gaussian kernels. As a result, the ratio for k- 
sharp kernel as well as Gaussian is R{ — 1.4i?vv. Ratio for the 
constraining scale is assumed the same, accordingly to the kernel 
used. The constraining mass is defined similarly to the filtered one: 
Mh = (47r/3) Qn^PoRi- 

For the single constraint with H — Wth we can see that the 
value of the constraining functional is bound to the mass associated 
with the host perturbation and to the averaged overdensity: 



A(7?h) 



C 



3Mh 



AnQn^poRl 



1 



(39) 



(note the remark on constraining scale made in the previous para- 
graph). This quantity carries the growth factor of the host perturba- 
tion for a given epoch, as a multiplier. 

Choice of the constraining kernel as well as filtering one is 
not unique and motivated by convenience usually. Top-hat kernel 
is a natural choice for a mass constraint as it has no 'tails' in the 
real space. On the oth er hand the k-sha rp filter suits better for ex- 
cursion set formalism (lBondetal.ll99l h. Also k-sharp kernel gives 
simple expressions in momentum space. Indeed, consider the sin- 
gle constraint with H = W = Wks- Denoting S*-"^ an uncon- 
strained variance (Eq. |6) we get the trivial relation for subhaloes 
(-Rw < -Rh): 



S = S^°'>{Rw)-S^">{Rh 



(0)/ 



(40) 



Something similar should take place for other kernels, although in 
a case of a larger number of constraints it might not be so. At Fig.[T] 
is presented the function Shim) for different kernels and different 
numbers of constraints. The models were computed for the host 
halo mass Mh = 10^^/i~^Mq. As we see the Top-hat kernel lead 
to the variance stable with respect to adding more constraints. Also 
the k-sharp has the same property for masses m < Mh but it gives 
oscillations for greater masses, seemingly due to oscillations in the 
real space. In the case of the Gaussian filter we see a more logical 
behaviour: the increase in the number of statistical constraints lead 
to the decrease in the variance in total. Note, since the function 
Sh ijn) is monotonically decreases, the lesser values of the variance 
cause the lesser amount of sub-haloes of associated mass, shifting 
mass PDF to the lower masses while the total amount of mass keeps 
unchanged and equals to Mh. Nevertheless, the adopting of Top- 
hat kernel instead of k-sharp, or Gaussian instead of Top-hat lead 
to effectively lower host mass and to younger sub-halo population. 
Applying different kernels for constraint and filter is probably not 
a good idea since kernels may interfere giving oscillations in real 
space. 



2.4 Linear growing of perturbations 

A linear evolution law for density perturbations i s known for an 
arbitr ary uniform and isotropic cosmology model (iBildhauer et al.l 
I1992I) . It depends on global density parameters of the Universe and 
gives the amplitude of decaying or growing modes as a function of 
scale factor or redshift. The solution for a growing mode is 



where 



-3/2 



X(a) 



db 



X(a) = ( 1 + 7^"+ 



1/2 



(41) 



(42) 
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Three constraints 



10* 10^' 10** 10"' 10^2 10" li)^^ 

m [/i^^Mg] 

Figure 1. Linearly evolved variance calculated for vaiious number of con- 
straints, from top to bottom: one (just kernel), two (the kernel and its 2- 
nd derivation) and three (the kemel, its 2-nd & 4-th derivations). Host halo 
mass is lO*^/i~^M0. Kernels used are: k-sharp (solid line). Top-hat (long- 
dashed) and Gaussian (short-dashed). Upper dotted line is the unconstrained 
variance, and horizontal dotted line mai'ks the level of nonlinearity, S\j = \. 
As seen, the Top-hat kemel gives variance which is invariant to the change 
of the number of constraints. For k-sharp kemel this property is tiue only for 
masses lesser than host. In both cases, the nonlinear mass is of order of the 
host halo mass. Using Gaussian kemel leads to the substantial dependence 
of the variance on the constraints number 

The evolution of perturbations inside an overdense or under- 
dense region can be examined in a pseudo-cosmological notation, 
if par ameters of sucii notation are cliosen appropriately dPeeblesI 
Il993l) . Namely, the need to set the density of the matter and dark 
energy to corresponding values in a host halo or void, then do the 
integral for growth factor using the new density parameters : 

= (1 + A)a^ 

= (43) 

1 — — 

Growth factor calculated in thise case we will denote D. Here A 
is the overdensity of a host perturbation linearly evolved to the 
present day, i.e. using the notation of Eq. |9l it can be written as 
D'*^-* = 0) Al. It is natural to bind this value to the overdensity 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



a 

Figure 2. Linear growth factor as the function of scale factor Perturbations 
are settled on a background with different overdensity A. The case A = 
corresponds to the unperturbed background. The solid line is the upper limit 
of the dependence when D(z = 0) = 5 (see text). 

averaged over the volume of the host perturbation, Eq. |39] namely 
A = A(_Rh), understanding the last as evolved to the present time. 
Substituting these definitions into equations above we get a law of 
growth on the given background. Properties of solutions are fol- 
lowing. If the linear overdensity of the host is positive, then the 
amplitude will grow faster turning to a greater value at present 
time. However, an arbitrarily large value of A can lead to only fi- 
nite linear amplitude a t present time with limiting value i3(0) = 5 
(Bildhauer et alJlT992L eq. 21 for f2m 2> 1). Otherwise, the arbi- 
trarily small value of the linear overdensity lead to smaller values 
of amplitude at present time, giving -D(O) = as a limiting case 
for A = -1. 



3 APPLICATIONS 

Our theory involves several parameters in addition to those in the 
original Press-Schechter approach. The parameters are the spatial 
scale of constraints i?h or associated mass A/h, and the mean over- 
density A. Constraining kemel can also be chosen fairly free. As 
we seen in Subsec. 12.31 the, different kernels with the same spa- 
tial scale effectively correspond to the different masses of the host 
halo, since they give different rates of decrease of the variance with 
spatial scale. Such an effective mass or scale can be used to cali- 
brate i?h for various kernels. The effective value of this parameter 
can be calculated as the spatial momentum of A''^'(r) profile or 
so. Other way is to bind the effective scale to the scale where con- 
strained variance crosses unity. Seems reasonable to use the model 
with H = W — Wth as a reference because the relations of the 
kind Eq.[39]looks naturally in the real space for this case. In the fol- 
lowing examples of applications we'll use the Top-hat kemel both 
for filter and for constraints. 

We need to check the correctness of our theory for various 
types of the real objects. As was mentioned in Subsec. 12.11 the 
stmcture of the host halo may affect the mass function of the sub- 
halo population. This impact should take place through the mass 
segregation and to be reflected in the sub-halo merging process. 
It is reliably enough to assume that such situation occurs in viri- 
alized haloes only. Hence, the knowingly non-virialized stmctures 
such as voids a nd friable supercluste rs, having low overdensities 
(A < 100, see lLacev&Coi3 ( ll993h ). are free from this effect of 
the 'stmcture biasing'. Last means that for these objects the ujm 
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factor in Eq.[TO]equals to unity and number density PDF becomes 



dm 



fn 



dS 



dm 



fs , 



(44) 



where fs from Eq.|2] The cumulative distribution function is 



7V(> m) 



dm — /™ 



(45) 



where is the total mass of an area of interest, i.e. of the 
sub-area of a supercluster or void. This value can be bound to 
the observed scale Rs at the given redshift via relation Ms = 
(47r/3) ilniPo^s (1 + Without resorting to calculations we can 
find qualitative properties of these distributions. Since the growth 
factor is limited for all possible values of the host overdensity, and 
the variance depends on the mass continuously (if constraining ker- 
nels are continuous), there are the limiting distribution finction for 
given constraining scale i?h or 7\'/h, corresponding to the limit of 
the evolution of sub-haloes population when the high-mass end of 
the PDF is the heaviest possible. Other limiting property is the low- 
mass end of the distribution where PDF tends to the unconstrained 
law which is rescaled accordingly to the growth factor value and 
independent of i\fh. 

Below we will very briefly touch three possible applications of 
the theory: the mass function in superclusters and also in the voids, 
and the bias relation for haloes in non-virialized hosts. The primary 
goal here to show the dependence of the distribution functions on 
the host overdensity. 



3.1 Superclusters 

Superclusters exhibit a very broad distribution of sizes from lQh~^ 
Mpc extending up to 150/i"^ Mpc, and mass up to lO^^/i^^M©. 
These objects form mainly in shells and filaments, although lumps 
are also observed. Clusters of galaxies are the characteristic compo- 
nents of superclusters. The typical dynamical mass of rich clusters 
is about 10^'' — lQ^^h~^yiQ. A number of rich clusters in super- 
cluster can vary from se veral to over t en, and their mass fraction 
can be about 50 per cent ( lBahcalllll999h . 

Superclusters can be characterized with density enhancement 
factor which is the number of objects (clusters or galaxies) in the 
super cluster related to the mean number of objects in the same vol- 
ume teahcall & Soneirairi984) . On the theoretical base the factor 
can be represented as the relation 



jV(>m) 



(46) 



where m is the lowest mass of objects we are counting. Here the Eq. 
I45l is used with Als = Afh, and the factor 1 + A is applied to trans- 
foiTn the Lagrange volume of the evolved host perturbation to the 
unperturbed Eulerian one. The density enhancement factor keeps 
information about total mass or radius and depth of the host per- 
turbation. On the Fig.[3]this function is presented for Afh — 10^'^ 
andlO^'';i"^M0,andfor various overdensity values. For the case 
if the host mass much exceeds the nonlinear mass (i.e. when the 
constraints role are negligible) the low-mass limit of last equation 
(m < lO^'^ft^^M©, which corresponds to the dwarf galaxies) al- 
lows to reduce the density enhancement factor for a fixed epoch to 
dependency on A only. 



E : 



— (1 + A) 




W 10" 10" 



101 



101 



Figure 3. Density enhancement factor as the function of the lowest sub-halo 
mass used for counting, and the host overdensity. Top panel: Mg = Mji = 
IQlSfe-lM©. Bottom panel: = = W^^H-^Mq. As seen, there 
is no dependense of the factor on the host halo mass or on the lowest sub- 
halo mass for small values of last, although the dependense is strong for 
moderate and high masses. In theory with unconstrained statistics we have 
E = 1 (not shown). 



so measuring E by counting objects of small masses we can re- 
store profile of the host overdensity of largest superclusters. How- 
ever, this may be difficult since the number of dwarf galaxies in 
a superscuster can be tens and hundreds of thousands, and it may 
be biased by a sub-clustering. In order to obtain the density profile 
by counting objects of moderate masses the following way might 
be used. Considering an external radius of a supercluster _Rs, count 
the number and the density enhancement factor of objects having 
mass greater than some value. Then solve Eo s 145 1 and 146 1 as sum- 
ing Ms corresponding to i?s (see the prevoius Subsec). Solving 
equations we get the mean overdensity A and the lowest mass used 
previously for counting, i.e. it is possible to calibrate the assumed 
lowest mass. The parameters found can be used further for smaller 
radii. 

Example of cumulative mass functions for the host Mh = 
W^h-^MQ is shown at Fig.gl We see that the mass of the most 
massive sub-halo (level A'^ = 1) increases several times when over- 
density changes from 3 to 20 but does not become higher than ~ 20 
per cent of the mass of host. On these figures and on the graph of 
the cumulative mass fraction of sub-population (Fig.|5} is seen how 
the distribution of sub-halo masses slides toward the high mass end 
when increasing host overdensity. This bias will be examined more 
detailed in Subsec. [33] 



(47) 



3.2 Voids 

Voids are formed from an underdensed regions by ousting the mat- 
ter out and building the cosmological walls herewith. Voids are 
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IQi" IQii IQi^ 1013 10" IQi^ lOi'' 

m [/i-^Mq] 

Figure 4. Cumulative distribution functions for sub-haloes inside a host. 
Top panel: Ms = = IO^^H-'^Mq. Bottom panel: Mg = Mh = 
IO^^/i^^Mq. The increment of the overdensity induces the redistribution 
of the mass function to higher mass values. Inset shows the same graphs 
zoomed in the high-mass end. The mass function for unconstrained case is 
very close to the case A = 0, as expected. 



1 




Figure 5. Cumulative distribution function of the mass fraction of sub- 
haloes inside a host Ms = Mjj = IO^^/i^^Mq for various linear over- 
densities. Designations of overdensities by lines are the same as in Fig.|4] 

charac terized by the typical s ize of tens megaparsecs and prolate 
shape jposter & Nelsonll2009l) . According to simulations, the over- 
densi ty of the matter inside a void can be as lo w as —0.95 at Hubble 
time jvan de Weygaert & van Kampen| [l993h . 

IPustilnik & Tepli^ovai 1201 10 showed that galaxies in the 
Lynx-Cancer void can form groups and filaments. So the sub-halo 
population in voids is not spatially smooth but undergoes cluster- 
ing. The distribution functions of sub-halo population inside voids 
are represented on Figs.|4]|5]and|6]for A = —0.75. 

Despite the fact that cosmological voids don't virialize they 




IQio 10^1 IQi" 10^=^ 10" 10^" 1016 
m [/i-iMg] 



Figure 6. Distribution function of the number density of sub-haloes in- 
side a host. Top panel: Mi^ = 10^^h~^yiQ. Bottom panel: Afjj = 
10^^'h~^M.Q. Designations of overdensities by lines are the same as in 
Fig.E] 



become nonlinear and the void's boundaries become walls, so the 
total amount of mass occupied by the structure distinguished as 
a void is by few orders smaller than mass of the host perturba- 
tion. As the mean overdensity in a void is negative, the sub-halo 
population should be younger comparing to the supercluster or to 
the unconstrained case, i.e. the mass function should be heavier in 
the low-mass end. For these reasons the deficit of massive galax- 
ies should be observed leading to sharper decay of the luminos- 
ity function for galaxies of moderate and high masses. This may 
be seen clearer on the number density distribution function, Fig. 
|6] Comparing to the original Press-Schechter theory without con- 
straints (dotted line) the PDF inside the void (long-dash-dotted line) 
has more sharp dec line, so the differents become significant for 
m > lOi^fe'iMfn. lPustilnik^ TepUakova (2011) didn't observe 
this effect for their Lynx-Cancer void galaxies sample. This can 
be explained by the fact that the sample is limited by the abso- 
lute magnitude ~ — 18 which correspon ds to dwarf galaxies hav- 
ing photometric mass M25 ~ 10^°Mq jKarachentsev et alj|2004 
the Hubble parameter is ft = 0.72), so the halo mass doesn't ex- 
ceed 10"Mq. This mass limit reveals the trivial behaviour of the 
low-mass end of mass function (Fig.[6]l. To distinguish the features 
our theory brings, the galaxies with absolute magnitude lesser than 
—20 or AI25 ^ 10" Mq should be observed. For example, if the 
host perturbation has the mass of order A/h ~ lO^ft^^Me and the 
mean overdensity of the void is A = —0.75 inside the sphere Ra = 
15h~^ Mpc (see simulations of Ivan de Weygaert & van Kampenl 
( Il993l) ). then the enclosed mass is Ms ~ lO^/i'^Mo. Such void 
should contain a few haloes of mass ~ 
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3.3 Bias factor 

Galaxy biasing is known effect when the bright gala xies tend to 
form near the high density peaks dBardeen etal.lll986h . In the dark 
matter section this property reflects the fact that density peaks 
are more clustered than underlying mass dMo & W hite 1996). In 
other words, the dark haloes trace the features of the background 
density field. Halo biasing was examined theoretically first by 
lBardeenetaIl l ll986t) . O n base of the Press- Schechter theory the 
bias was obtained first bv lMo & Whit3 l ll99 ^ as a relation between 
correlation f unctions of haloes and m ass. Their ex pression was then 
extended by ISheth & TormenI ( Il999l) . revised by Ijind ( Il998h and 
|Porcianietal.l ( ll999h . as well adjusted with numerical N-body re- 
sults. 



In approach of ISheth & TormenI ( ll999 S the Eulerian bias fac- 
tor was defined as the linear part of the relative overabundance of 
the halo number inside volume V containing given mass M from 
the mean halo number in the same volume: 

N(ni\M, V) 



3hh 



1 



(48) 



n(o)(m)F 

Present day value of the factor obtained in the Press-Schechter the- 
ory is 

5, D("'> 



= 1 + 



(49) 



Ijin^d 19981) and lPorcianietal .Id 19991) faced the good agreement be- 
tween this function and accurate numerical measurements for large 
halo masses, but found a departure when the mass of the halo is 
small; the model of Mo and White underestimated biasing for lower 
mass limit leading to a constant value of the bias factor while nu- 
merical measurements revealed a rising of it with mass decreasing. 

The bias factor of a kind Eg. 149 1 depends on mass of the halo 
but not on the background. In the previous numerical works the 
procedure of calculating bias was constructed in such a way that 
the factor was determined being averaged on a big enough compu- 
tational volume. Using formulation of the mass function proposed 
in this paper we can estimate the biasing effect stipulated by only 
the host halo mass and overdensity. In a manner Eq.|48]the relative 
overabundance of the halo number can be written as 



(50) 



where the volume tranformation factor 1 + A is considered, and 
'zero' index means the analogous values from the unconstrained 
theory. Eulerian bias factor then defined as 

, _ d5hh 



dA 



(51) 



In this foiTnulation both the dynamical age of the sub-halo popula- 
tion (determined by host overdensity A) and power spectrum con- 
straints (determined by host mass Mh) affect the bias. On the Fig.|7] 
presented the bias factor for Mh = lO^^/i~^M0 and various host 
overdensities. The more rarefied background hosts younger sub- 
halo population which becomes more biased in a lower masses. At 
the same time the dense background demonstrates very strong bi- 
asing in a high masses comparing to the approximating formula 
of Mo & White. This discrepancy may be cleaned by the follow- 
ing way. The large scale structures of different overdensities have 
the different volume filling factor in the Universe. Namely, voids 
having negative overdensity occupy the major part of the Universe 
while the walls, filaments and other structures having positive over- 
density occupy the volume with factor of order few per cents. Con- 
sidering this it may be possible to match the rising low-mass limit 
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Figure 7. Present day Eulerian bias factor for My^ = IO^^'/i^^Mq and 
various host overdensities. It is seen, the younger sub-halo population (A is 
lower), the more it is bia sed for the lower ma sses of sub-haloes. Dotted line 
shows the bias factor bv lMo&Wh"it3il996l) . 



of the bias djing|ll99^ while to attenuate an exponential growth of 
the bias factor in the high-mass limit. 



4 DISCUSSION 

Comparing with the original Press-Schechter theory, our theoreti- 
cal model has additional parameters: mean overdensity of the host 
perturbation A, its spatial scale iih, constraining kernel H defin- 
ing profile of the host perturbation, and its spatial momenta Ca, 
first of which is A. From the original theory, our approach is for- 
mally characterized in using of different mass-variance function 
and assuming some different background cosmology. Thereby, all 
the statistics it gives is spatially uniform. In addition, neither host 
perturbation profile nor its overdensity changes the merging pro- 
cess of the sub-halo population. That is, our theory in its present 
form is not able to model properly the virialized haloes, e.g. haloes 
of galaxies and galaxy clusters. As was mentioned above, the virial- 
ization process may result in mass segregation which alternates the 
merger rates for low- and-high-mass sub-haloes. The process has 
more significant influence on the high-mass end of the distribution 
function. This effect can be named the 'structure biasing'. 

Correction to the distribution functions accounting to the 
structure biasing may have form of the multiplier, like in Eg. 1101 
depending on parameters noted above, and also on sub-halo mass. 
It is possible to estimate the high-mass limit of this multiplier using 
virial considerations. In general case the spatial distribution of the 
density per unit mass interval is 

dpr, 
dm 



(52) 



where uj = aj(r, m) is a factor responding for the structure bias. 
Cumulative mass of a sub-halo population inside a sphere r is 



M(< r) = An / dr' r 



dm ^ 

am 



= 3Mh 



da; X 



(53) 



dmujfn 



where x — r/ Rh is introduced. The lj factor depends on the viri- 
alization process and reflects as density profile of a host halo so the 
sub-haloes mass segregation inside it. Without going into details 
of virialization we can denote that sub-haloes of high mass locate 
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naturally near the centre while low-mass haloes occupy the entire 
host volume. The meaning of ljJ then might be the relative volume 
filled by sub-haloes of a given mass. The factor should decrease 
with sub-halo mass increasing, from a constant value of order unity 
or greater, for m <C Mh, to a somewhat small value in the limit 
m — > Afh. The detailed structure of the high mass limit can be ex- 
amined as follows. F or definiteness let th e virialized host halo has 
NFW density profile jNavarro et al.ll997b : 

Acpo 



Pnfw 



(54) 



(1 + r/vs)^ r/ra 

where Vs is the scale radius; Ac is the characteristic overdensity. 
The mass profile with it's Taylor expansion is 



Afh 



3Ac 

200c3 V 2r| 



(55) 



where c is the concentration factor, and i\fh has the meaning 
of the mass inside a sphere with overdensity 200, i.e. A/h = 
(47r/3) 200^0^200. ^^m-Rg = 200r^oo. and r2oo = cr^. In the 
vicinity of the centre the heavy mass sub-haloes dominate by mass 
so we can assume mass profile of the population (Eg. 153) and of the 
host halo (Eq.|55t matching there. Hence, the structure bias factor 
can be expressed in the form 

f^m r-/rs 

with the natural mass conservation claim 

drriLJmfm = 1 • (57) 

The mass part Um of the structure bias factor can be obtained using 
virial theorem and equipartition. Gravitational potential of NFW 
halo is 



3Ac GA/h ln(l + r/rs) 



200c3 



r/r 



(58) 



Write the virial theorem, neglecting an angular momentum of the 
motion, then expand the virial at the lowest radii: 



dr 



+ 



(59) 



Assuming equipartition at the lowest energies, we get (v^) oc m ^, 
and hence (r) oc . Now define CJ™ as the volume factor, LOm oc 
(r)^, then 



ijJm OC — - for m ~ Mh 



(60) 



As we see, the structure bias factor reduces the amount of massive 
sub-haloes simultaneously increasing the amount of small ones. 



5 CONCLUSION 

We proposed the modification of the Press-Schechter theory of 
dark halo clustering. The basis of modification is the use of the 
statistical constraints imposed on the initial gaussian random fluc- 
tuations. Constraints has the form of the linear functionals of the 
fluctuations' overdensity field. Kernels of constraints determine the 
shape and, mainly, the spatial scale of the host perturbation which 
evolves to the supercluster or void. The resulting non-gaussianity 
of the random field leads to existing of the upper limit of mass of 



the sub-haloes formed inside the host perturbation. Amplitude of 
the perturbation also determines the age of the sub-halo population 
in the sense of the merging process. 

Since the only spherical host perturbation is adopted in this 
paper we can not apply our theory to the variety of morphologies 
of superclusters but to the lumps only. However, it is not difficult to 
generalize the approach used here to an arbitrary shape of the host 
perturbation to the wall or filament. 

Calculations show, the distribution function for sub-haloes 
with mass fraction 10^* — 10^'' of the host and greater reveal de- 
viations from the original Press-Schechter predictions, while devi- 
ations are negligible for lower masses so they can be expressed as a 
rescaling. If the host halo mass is A4 = lO^^/i"^M0 audits over- 
density is A = 20, the deviations can be detected for sub-haloes 
of the mass W^^ H'^Mq and greater. Assuming dark-to-light mass 
relation be 10, such haloes correspo nd to the galaxies having ab- 
solute magnitude —20 and brighter jKarachentsev et al]|2004h . To 
check the theory in a voids' galaxies is necessary to count ones with 
absolute magnitudes lesser than —20. However, there may be only 
a few such galaxies inside a tens-megaparsecs void. Also using the- 
ory, it is possible to recall the density profile of the superclusters. 

Approach used in this paper does not allow to predict the dis- 
tribution functions in virialized haloes. However, in previous sec- 
tion showed, the effect of the structure biasing should compen- 
sate the shift of the distribution function inside an overdense halo 
toward the high masses. This can explain the good fitting of the 
Press-Schechter function to a wide range of observations and nu- 
merical experiments. 
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